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The low temperature unitary Bose gas is a fundamental paradigm in few-body and many-body 
physics, attracting wide theoretical and experimental interest. Here we first present a theoretical 
model that describes the dynamic competition between two-body evaporation and three-body re¬ 
combination in a harmonically trapped unitary atomic gas above the condensation temperature. We 
identify a universal magic trap depth where, within some parameter range, evaporative cooling is 
balanced by recombination heating and the gas temperature stays constant. Our model is developed 
for the usual three-dimensional evaporation regime as well as the 2D evaporation case. Experiments 
performed with unitary ^®^Cs and ^Li atoms fully support our predictions and enable quantitative 
measurements of the 3-body recombination rate in the low temperature domain. In particular, we 
measure for the first time the Efimov inelasticity parameter rj, — 0.098(7) for the 47.8-G d-wave 
Feshbach resonance in ^^®Cs. Combined '^^^Cs and ^Li experimental data allow investigations of loss 
dynamics over two orders of magnitude in temperature and four orders of magnitude in three-body 
loss. We confirm the 1/T^ temperature universality law up to the constant rjt. 

PACS numbers: 05.30.Jp Boson systems 

05.70.Ln Nonequilibrium and irreversible thermodynamics 

34.50.-s Scattering of atoms and molecules 

51.30.-|-i Thermodynamic properties, equations of state 


I. INTRODUCTION 

Resonantly interacting Bose systems realized in ultra¬ 
cold atomic gases are attracting growing attention thanks 
to being among the most fundamental systems in na¬ 
ture and also among the least studied. Recent theoret¬ 
ical studies have included hypothetical BEC-BCS type 
transitions [1-5] and, at unitarity, calculations of the uni¬ 
versal constant connecting the total energy of the system 
with the only energy scale left when the scattering length 
diverges: jra [6-9]. The latter assumption 

itself remains a hypothesis as the Efimov effect might 
break the continuous scaling invariance of the unitary 
Bose gas and introduce another relevant energy scale to 
the problem. A rich phase diagram of the hypothetical 
unitary Bose gas at finite temperature has also been pre¬ 
dicted [10, 11]. 

In experiments, several advances in the study of the 
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resonantly interacting Bose gas have recently been made 
using the tunability of the s-wave scattering length a near 
a Feshbach resonance. The JILA group showed signa¬ 
tures of beyond-mean-field effects in two-photon Bragg 
spectroscopy performed on a ®®Rb BEG [12], and the 
ENS group quantitatively studied the beyond mean-field 
Lee-Huang-Yang corrections to the ground state energy 
of the Bose-Einstein condensate [13]. Logarithmic be¬ 
havior of a strongly interacting 2D superfluid was also 
reported by the Ghicago group [14]. Experiments have 
also started to probe the regime of unitarity (1/a = 0 
directly. Three-body recombination rates in the non¬ 
degenerate regime have been measured in two different 
species, ^Li [15] and [16], and clarified the temper¬ 
ature dependence of the unitary Bose gas lifetime. In 
another experiment, fast and non-adiabatic projection of 
the BEG on the regime of unitarity revealed the establish¬ 
ment of thermal quasi-equilibrium on a time scale faster 
than inelastic losses [17]. 

In a three-body recombination process three atoms col¬ 
lide and form a dimer, the binding energy of which is 
transferred into kinetic energies of the colliding partners. 
The binding energy is usually larger than the trap depth 
and thus leads to the loss of all three atoms. Because 
three-body recombination occurs more frequently at the 
center of the trap, this process is associated with “anti- 
evaporative” heating (loss of atoms with small poten¬ 
tial energy) which competes with two-body evaporation 
and leads to a non trivial time dependence for the sam¬ 
ple temperature. In this paper, we develop a theoretical 
model that describes these atom loss dynamics. We si- 
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multaneously take into account two and three-body losses 
to quantitatively determine each of these contributions. 
We predict the existence of a magic value for the trap- 
depth-over-temperature ratio where residual evaporation 
compensates for three-body loss heating and maintains 
the gas temperature constant within some range of pa¬ 
rameters. We then apply our model to analyze the loss 
dynamics of ^^^Cs and ^Li unitary Bose gases prepared 
at various temperatures and atom numbers. Compar¬ 
ing measurements in these two different atomic species 
we find the dynamics to be universal, i.e. in both sys¬ 
tems the three-body loss rate is found to scale universally 
with temperature. Excellent agreement between theory 
and experiment confirms that the dynamic evolution of 
the unitary Bose gas above the condensation tempera¬ 
ture can be well modelled by the combination of two and 
three-body interaction processes. 


II. MODEL 

A former study developed for measuring three-body 
decay in trapped ^^^Cs [18] atoms has proposed a model 
to describe the time evolution of the atom number N and 
the temperature T taking into account the three-body re¬ 
combination induced loss and the heating associated with 
it. This model is valid in the limit of deep trapping poten¬ 
tials (trapping depth much larger than the atom’s tem¬ 
perature) and for temperature independent losses. Here 
we generalize this model to include evaporation induced 
cooling and the associated atom loss, as well as the tem¬ 
perature dependence of the three-body loss rate. 


A. Rate equation for atom number 

The locally defined three-body recombination rate 
L^n^(r)/?> leads, through integration over the whole vol¬ 
ume, to the loss rate of atoms: 


where fcth = y/mkBT / h, Rq is the three-body parame¬ 
ter, and the Efimov inelasticity parameter 77 * character¬ 
izes the strength of the short range inelastic processes. 
Here, h is the reduced Planck’s constant, is the Boltz¬ 
mann’s constant, and sq = 1.00624 for three identical 
bosons [19]. The matrix element sn relates the incom¬ 
ing to outgoing wave amplitudes in the Efimov scattering 
channel and shows the emerging discrete scaling symme¬ 
try in the problem (see for example Ref. [20]). Details are 
given in the supplementary material to Ref. [15] for the 
calculation of Sii(fca), where a is the scattering length 
and k is the relative wavenumber of the colliding part¬ 
ners. Because of its numerically small value for three 
identical bosons at unitarity, we can set |sii| = 0 and L 3 
is well approximated by: 


L 3 
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{kBTy 
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where A 3 is a temperature-independent constant. Assum¬ 
ing a harmonic trapping potential, we directly express 
the average square density (n^) through N and T. In 
combination with Eq. (3), Eq. (1) is represented as: 


where 
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with w being the geometric mean of the angular frequen¬ 
cies in the trap. 

To model the loss of atoms induced by evaporation, we 
consider time evolution of the phase-space density distri¬ 
bution of a classical gas: 


-U(r)/kBT -p'^/2mkBT 

" (27rA)3 


( 6 ) 


dN 

di 
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dV = -L3{n^)N, 


( 1 ) 


where the factor of 3 in front on the integral reflects the 
fact that all 3 atoms are lost per each recombination 
event. In the following, we neglect single-atom losses 
due to collisions with the background gas and we assume 
that two-body inelastic collisions are forbidden, a condi¬ 
tion which is fulfilled for atoms polarized in the absolute 
ground state. 

An expression for the three-body recombination loss 
coefficient at unitarity for a non-degenerate gas has been 
developed in Ref. [15]. Averaged over the thermal distri¬ 
bution it reads: 


L 3 — 


72y37r2A(l-e-4''*) 




(l—|sii|^)e 

|1 -b {kRo)-^'^^°e-‘^-^‘Sii\'^ 


( 2 ) 


which obeys the Boltzmann equation. Here uq is the 
central peak density of atoms, AdB = (27r/i^/mfcBT)^^^ 
is the thermal de Broglie wavelength, and U{r) is the 
external trapping potential. The normalization con¬ 
stant is fixed by the total number of atoms, such that 
/ /(r, p)d3pd3r = N. 

If the gas is trapped in a 3-D trap with a potential 
depth U, the collision integral in the Boltzmann equation 
can be evaluated analytically [21]. Indeed, the low-energy 
collisional cross-section 

reduces at unitarity to a simple dependence on the rel¬ 
ative momentum of colliding partners: crik) = Sir/k'^. 
However, not every collision leads to a loss of atoms due 
to evaporation. Consider 

77 = U/kBT. 


X 
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( 8 ) 
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In the case of t] » 1, such loss is associated with a 
transfer of large amount of energy to the atom which 
ultimately leads to the energy independent cross-section. 
This can be understood with a simple argument [22] . As¬ 
sume that two atoms collide with the initial momenta pi 
and P 2 - After the collision they emerge with the mo¬ 
menta ps and P 4 , and if one of them acquires a mo¬ 
mentum jpsj > \/2mU. Then, |p 4 | is necessarily smaller 
than the most probable momentum of atoms in the gas 
and jpaj |p 4 |. In the center of mass coordinates the 
absolute value of the relative momentum is preserved, 
so that 5 IP 1 -P 2 I = l\P3 - P4\ « sIpsI- Assuming 
|p 3 1 = y/‘2mU, we get |pi — p 2 1 = V^mU. Substituting 
the relative momentum in the center of mass coordinate, 
hk = ^{pi — P 2 ), to the unitary form of the collisional 
cross-section, we find the latter is energy independent: 


O'er 


167rh^ 

mU ’ 


( 9 ) 


and the rate-equation for the atom number can be writ¬ 
ten as: 


=-TeviV, rev = no(T(7ue (10) 

at Ve 

The peak density is no = N/Ve, where I 4 is the effec¬ 
tive volume of the sample. In the harmonic trap V). can 

be related to uJ and the temperature T: 14 = ■ 

The ratio of the evaporative and effective volumes is de¬ 
fined by [21]: 

^ = »7 - 4i? (3,77), (11) 

where R{a,r]) = and P{a,r]) is the incomplete 

Gamma function 


P{a,v) 




Finally, taking into account both three-body recombi¬ 
nation loss (see Eqs. (4),(5)) and evaporative loss, we can 
express the total atom number loss rate equation as: 


dA^ _„KviV2 

dt ~ TS 14 T ’ 


( 12 ) 


where 


16 

TT k^U 


(13) 


Note that 77 and the ratio of the evaporative and effec¬ 
tive volumes explicitly depend on temperature and 72 is 
temperature independent. 


B. Rate equation for temperature 

1. ‘Anti-evaporation’ and recombination heating 

Ref. [18] points out that in each three-body recom¬ 
bination event a loss of an atom is associated with an 
excess of k^^T of energy that remains in the sample. 
This mechanism is caused by the fact that recombina¬ 
tion events occur mainly at the center of the trap where 
the density of atoms is highest and it is known as ‘anti¬ 
evaporation’ heating. We now show that the unitary 
limit is more ‘anti-evaporative’ than the regime of finite 
scattering lengths considered in ref. [18] where L 3 is tem¬ 
perature independent. We separate center of mass and 
relative motions of the colliding atoms and express the 
total loss of energy per three-body recombination event 
as following: 


^3b = -J ((Ac„4 + 3t/(r)) 

+ ^(A 3 (fc)i?fe)|dV (14) 

The first two terms in the parenthesis represent the mean 
center-of-mass kinetic energy (Acm) = (A’cm)/ 2 A 1 and 
the local potential energy 3C/(r) per each recombination 
triple. M = 3m is the total mass of the three-body 
system. The last term stands for thermal averaging of 
the three-body coefficient over the relative kinetic energy 
Ek = {HkY/2p where p is the reduced mass. 

Averaging the kinetic energy of the center of mass mo¬ 
tion over the phase space density distribution (Eq. ( 6 )) 
gives (Acm) = I^bE. Then the integration over this term 
is straightforward and using Eq. (1) we have: 

- I ^^^^^{E,m)d^r = ^kBTN (15) 

The integration over the second term can be easily 
evaluated as well: 

/ or 1 

-^n^{r)U{r)d^r = -k^TN (16) 

To evaluate the third term we recall the averaged over 
the thermal distribution expression of the three-body re¬ 
combination rate in Eq. (2). Now its integrand has to 
be supplemented with the loss of the relative kinetic en¬ 
ergy per recombination event E^. Keeping the limit of 
Eq. (3) this averaging can be easily evaluated to give 
{Lo{k)Ek) = LoksT. Einally, the last term in Eq. (14) 
gives: 

- J '^{Loik)Ek)d^r = ^kBTN (17) 

Einally, getting together all the terms, the lost energy per 
lost atom in a three-body recombination event becomes: 
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This expression shows that unitarity limit is more 
‘anti-evaporative’ than the regime of finite scattering 
length {k\a\ < 1). As the mean energy per atom in the 
harmonic trap is Sk^T, at unitarity each escaped atom 
leaves behind (3 — 4/3)fcB7’ = {^/i)k^T of the excess 
energy as compared to Ik^T when L 3 is energy indepen¬ 
dent. In the latter case, thermal averaging of the relative 
kinetic energy gives {Ek) = ik'QT, thus Esh/N = 2k'QT. 

Eq. (18) is readily transformed into the rate equation 
for the rise of temperature per lost atom using the fact 
that Asb = iNk^T in the harmonic trap and Eq. (4): 


dT _ 5 T 
d^ “ 


(19) 


Another heating mechanism pointed out in Ref. [18] 
is associated with the creation of weakly bound dimers 
whose binding energy is smaller than the depth of the 
potential. In such a case, the three-body recombination 
products stay in the trap and the binding energy is con¬ 
verted into heat. 

In the unitary limit, this mechanism causes no heat¬ 
ing. In fact in this regime, as shown in the supplementary 
material to Ref. [15], the atoms and dimers are in chem¬ 
ical equilibrium with each other, e.g. the rate of dimer 
formation is equal to the dissociation rate. We therefore 
exclude this mechanism from our considerations. 



N/N 

in 


FIG. 1: N-T phase space representation of ‘anti-evaporation’ 
heating and evaporative cooling dynamics for different values 
of the initial rjin parameter. The “magic” rjm satisfies the 
condition dT /dN = 0. For lower values of rjin the “magic” 
is not reached during the evolution of the gas.The figure is 
drawn for experimental parameters of ^^®Cs atoms presented 
in Sec. Iff. For these conditions 77 m coincides with rjin ~ 8. 


2. Evaporative cooling 

“Anti-evaporative” heating can be compensated by 
evaporative cooling. The energy loss per evaporated 
atom is expressed as: 

E = N{r] + k)kBT (20) 


where, as before, the temperature dependence remains in 

T], 

Finally, combining the two processes of recombination 
heating (Eq. (19)) and evaporative cooling (Eq. (24)) we 
get: 


dt 


T 

3 


/5 

(^373^-726 


-V 




{r] + k 



(25) 


where k in a harmonic trap is [ 21 ]: 

' , P{5,v) Ve 

P (3,7?) Rev 


( 21 ) 


with 0 < k < 1 . 

In a harmonic trap, the average energy per atom is 
Sk^T = Taking the derivative of this equation and 
combining it with Eq. (20) we get: 

f N 

“ TV ^ ~ 

From Eqs. (10) and (22), evaporative cooling is expressed 
as: 


3—=-Tev (t?- h k - 3)T, (23) 

dt 

This equation can be presented in a similar manner as in 
the previous section: 


dT 

dt 


= -726 {r) + k 


3)^^ 

’ T 3' 


(24) 


Eqs. (12) and (25) form a set of coupled rate equations 
that describe the atom loss dynamics. 


C. N-T dynamics and the “magic” rjni 

To study atom number and temperature dynamics we 
solve Eqs. (12) and (25) numerically for different initial 
values of 77 , referred to as from here on. As an illus¬ 
tration, 72 and 73 are evaluated based on parameters of 
the ^^^Cs experiment discussed in Sec. III. The system 
dynamics in iV — T phase space is represented in Fig. 1. 
All represented values of rji^ lead to a decrease in temper¬ 
ature for small atom numbers indicating that evaporative 
cooling always wins for asymptotic times where the atom 
density becomes small. This weakens the three-body re¬ 
combination event rate and effectively extinguishes the 
heating mechanism altogether. Large values of the initial 
rjin cause initial heating of the system which is followed 
by a flattening of the temperature dependence at a cer¬ 
tain atom number (grey dashed and dark yellow dotted 
lines) that defines the “magic” 77 m. In Fig. 1 the solid 
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FIG. 2: Universal plot rjm vs a = N{hu;/ksT)^{1 — 

(hue curve). The blue solid circles correspond to the results 
obtained for ^^®Cs in Fig. 3(a) with rj^ = 0.098. The red 
solid squares correspond to the ^Li data of Fig. 3(c) with 
rj* = 0 . 21 . 


green line represents the special case when rym = rjm- Ex¬ 
perimentally, r]{T, N) is tuned to satisfy this special case 
for a given initial temperature and atom number. As 
it is seen in Fig. 1, lower initial values of can never 
reach the necessary condition for rjm in their subsequent 
dynamics (orange dotted-dashed line). 

The value of ri^{T, N) is found by solving the equation 
dT/dN = 0, i.e. when T{N) becomes independent on 
the atom number up to the first order in N. From the 
general structure of this equation, we see that is solely 
function of the dimensionless parameter 

Up to a factor (1 —depends only on the phase- 
space density N^huj/ksT)^ of the cloud. We plot in Fig. 2 
the dependence of r]^ vs a. Since our approach is valid 
only in the non-quantum degenerate regime where the 
momentum distribution is a Gaussian, we restricted the 
plot to small (and experimentally relevant) values of a. 


D. 2D evaporation 

The above model was developed to explain 3D isother¬ 
mal evaporation in a harmonic trap and experiments with 
^^^Cs presented below correspond to this situation. Our 
model can also be extended to 2D isothermal evapora¬ 
tion, as realized in the ^Li gas studied in Ref. [15] and 
presented below. In this setup, the atoms were trapped in 
a combined trap consisting of optical confinement in the 
radial direction and magnetic confinement in the axial di¬ 
rection. Evaporation was performed by lowering the laser 


beam power which did not lower the axial (essentially in¬ 
finite) trap depth due to the magnetic confinement. Such 
a scenario realizes a 2D evaporation scheme. Here, we ex¬ 
plore the consequences of having 2D evaporation. In the 
experimental section we will show the validity of these 
results with the evolution of a unitary ^Li gas. 

Lower dimensional evaporation is, in general, less ef¬ 
ficient than its 3D counterpart. ID evaporation can be 
nearly totally solved analytically and it has been an in¬ 
tense subject of interest in the context of evaporative 
cooling of magnetically trapped hydrogen atoms [21, 23, 
24]. In contrast, analytically solving the 2D evaporation 
scheme is infeasible in practice. It also poses a rather 
difficult questions considering ergodicity of motion in the 
trap [25]. The only practical way to treat 2D evapora¬ 
tion is Monte Carlo simulations which were performed 
in Ref. [25] to describe evaporation of an atomic beam. 
However, as noted in Ref. [25], these simulations follow 
amazingly well a simple theoretical consideration which 
leaves the evaporation dynamics as in 3D but introduces 
an ’effective’ rj parameter to take into account its 2D 
character. 

The consideration is as following. In the 3D evapo¬ 
ration model, the cutting energy Cc is introduced in the 
Heaviside function that is multiplied with the classical 
phase-space distribution of Eq. ( 6 ) [21]. For the 2D 
scheme this Heaviside function is Y{€c — ej,), where Cc 
is the 2D truncation energy and e_L is the radial energy 
of atoms in the trap, the only direction in which atoms 
can escape. Now we simply add and subtract the axial 
energy of atoms in the trap and introduce an effective 3D 
truncation energy as following: 

E(ec ~ C-l) = T((cc -I- Cz) — (e_L + Cz)) = Y (e[:® — etot); 

(27) 

where etot is the total energy of atoms in the trap and 
the effective truncation energy is given e®® = tc + — 

+ where we replaced Cz, by its mean value k^T in a 
harmonic trap. The model then suggests that the evapo¬ 
ration dynamics follows the same functional form as the 
well established 3D model, but requires a modification of 
the evaporation parameter ( 8 ): 

7?®® = 77 + 1 , (28) 

Then, the experimentally provided 2D 77 should be com¬ 
pared with the theoretically found 3D 77 ®® reduced by 1 
(i.e. 77 ®® — 1 ). 

III. EXPERIMENTS 

In this section, we present experimental T{N) trajec¬ 
tories of unitary and ^Li gases, and show that their 

dynamics are given by the coupled Eqns. (12) and (25). 
The ^^^Cs Feshbach resonance at 47.8 Gauss and the 
^Li Feshbach resonance at 737.8 Gauss have very sim¬ 
ilar resonance strength parameter Sj-es = 0.67 and 0.80 
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respectively [26, 27]) and are in the intermediate cou¬ 
pling regime (neither in the broad nor narrow resonance 
regime). We first confirm the existence of a “magic” rj^ 
for unitarity-limited losses for both species, with either 
3D or 2D evaporation. Then we will use the unitarity- 
limited three-body loss and the theory presented here to 
determine the Efimov inelasticity parameter of the nar¬ 
row 47.8-G resonance in ^^^Cs which was not measured 
before. 


A. N -T fits 

We prepare the initial samples at Tin and A^in as de¬ 
scribed in the Supplementary Materials. We measure the 
atom number N{t) and the temperature T{t) from m- 
situ absorption images taken after a variable hold time t. 
In Fig. 3(a), we present typical results for the evolution 
T{N) of the atom number and temperature of the gases, 
and we furthermore treat the hold time t as a param¬ 
eter. We also plot the relative temperature T/Tjn as a 
function of the relative atom number N/Ni^ for the same 
data in Fig. 3(b), and for ’’Li in Fig. 3(c). We then per¬ 
form a coupled least-squares fit of the atom number and 
temperature trajectories, Eqs. (12) and (25), to the data. 
We note that with our independent knowledge of the ge¬ 
ometric mean of the trapping frequencies, w, the only 
free fit parameters apart from initial temperature and 
atom number are the trap depth U and the temperature- 
independent loss constant A 3 . The solid lines are the 
fits (see Supplementary Materials) to our theory model, 
which describe the experimental data well for a large va¬ 
riety of initial temperatures, atom numbers and relative 
trap depth. We are able to experimentally realize the full 
predicted behavior of rising, falling and constant-to-first- 
order temperatures. 


B. Magic rj 

The existence of maxima in the T — N plots confirms 
the existence of a “magic” relative trap depth 77 ^, where 
the first-order time derivative of the sample temperature 
vanishes. Using the knowledge of 77 * for both ’^^Cs and 
"Li, we can compare the observed values of 771,1 to the pre¬ 
diction of Fig. 2 (note that in the case of ’’Li, we plot 
that enters into the effective 3D evaporation model). We 
see that for both the 3D evaporation ’^^Cs data and 2D 
evaporation ’Li data, the agreement between experiment 
and theory is remarkable. 

Furthermore, in the Supplemental Materials we show 
that from the three-body loss coefficients and the evap¬ 
oration model, we can predict the trap depth, which is 
found in good agreement with the value deduced from 
the laser power, beam waist, and atom polarizability. 



Number of atoms N 




FIG. 3: (Color online) Evolution of the unitary gas 

in (a) absolute and (b) relative nnmbers. The solid lines are 
fits of the data using the theory presented here, and the htted 
initial relative trap depth 771,1 = Ujk^T^VL is given in the legend. 
Error bars are statistical. The condition for (dr/dA'’)|t=o is 
expected for ? 7 in ~ r/m ~ 8 . 2 , very close to the measured data 
for 77 in = 7.4 (green lines in a) and b)). (c) Evolution of the 
unitary ’Li gas. The solid lines are hts of the data using 
our 2D evaporation model, and the fitted initial relative trap 
depth rjin = U/k^Ti^ is given in the legend. Error bars are 
statistical. In 2D evaporation, r]i-n « r]^ = 77 m -f 1 = 8.5 is 
required to meet the (dr/dA)|t=o condition, and is found in 
excellent agreement with the measured value 8.5 (green line 
in c)), see text. 
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C. Universality of the three-body loss 

As the last application we now show the validity of the 
L 3 oc law for the tree-body loss of unitary ^Li and 
^^^Cs Bose gases. Because both species are situated at 
the extreme ends of the (stable) alkaline group, they have 
a large mass ratio of 133/7 = 19 and the temperature 
range is varied over two orders of magnitude from 0.1 |xK 
to 10 pK. We determine the three-body loss coefficients 
A 3 from fits to decay curves such as shown in Fig. 3. We 
present in Fig. 4 the results for the rate coefficient L 3 , 
which varies over approximately two orders of magnitude 
for both species. In order to emphasize universality, the 
loss data is plotted as a function of {m/mu)^T^, where 
toh is the hydrogen mass. In this representation, the 
unitary limit for any species collapses to a single universal 
line (dotted line in Fig. 4, cf. Eq. (3)). 



FIG. 4: (Color online) The magnitude of three-body loss rate 
at unitarity for ^Li (red) and '^^^Cs (blue) with the respective 
±1 standard deviation (shaded areas). On the horizontal 
axis, masses are scaled to the hydrogen atom mass mu- The 
dashed line represents the unitary limit (Eq. (3) with rjt, —>■ 
00 ). Solid lines are predictions of universal theory [15] with 
rjt — 0.21 for ^Li and 77 * = 0.098(7) for ^^^Cs, see text. The 
data confirms the universality of the I /3 oc T~^ law. 

For ^Li, we cover the 1-10 pK temperature range. 
We find for the temperature-independent loss coefficient 
A 3 = 3.0(3) X 10“^° cm®pK^s“^, very close to the uni¬ 
tary limit A™®'^ « 2.7 x 10“^° cm®pK^s“^. It is also close 
to the value A 3 = 2.5(3) x 10“^° cm®pK^s“^ found in 
[15] with a restricted set of data, and to the predicion 
from Eq. (2) with 77 * = 0.21 from [28] (red solid line in 
Fig. 4). We cannot measure 77 * here because the ’^Li data 
coincides with the unitary limit. 

Furthermore the quality of the ^^^Cs temperature and 
atom number data enables us to directly measure the 
previously unknown 77 * parameter of the 47.8-G Fesh- 
bach resonance. The standard technique for obtaining 
77 * is measuring the three-body loss rate L 3 (a, T —0) as 
a function of scattering length in the zero-temperature 
limit, and subsequent fitting of the resulting spectrum to 


universal theory. However, for a given experimental mag¬ 
netic field stability, this method becomes hard to put into 
practice for narrow resonances like the 47.8-G resonance 
in ^^^Gs. Instead, we use the fits to our theory model in 
order to obtain 77 * from A 3 . We cover the 0.1-1 pK range 
and find A 3 = 1.27(7) x 10“^^ cm®pK^s“^. Plugging this 
number into Eq. (3) , we deduce a value for the Efimov 
inelasticity parameter 77 * = 0.098(7). The corresponding 
curve is the blue line in Fig. 4 and is significantly be¬ 
low the unitary line because of the smallness of 77 *. This 
new value is comparable to the Efimov inelasticity pa¬ 
rameter found for other resonances in ^^^Gs, in the range 
0.06...0.19 [29, 30]. 

The plot of the full theoretical expression Eq. (2) for 
in Fig. 4 (full lines) requires an additional pa¬ 
rameter describing three-body scattering around this Fes- 
hbach resonance, the so-called three-body parameter. It 
can be represented by the location of the first Efimov 
resonance position [31]. Because of the lack of ex¬ 
perimental knowledge for the 47.8-G resonance, we take 
the quasi-universal value = —9.73(3)rvdWj ?’vdw be¬ 
ing the van-der-Waals radius, for which theoretical ex¬ 
planations have been given recently [31-33]. The theory 
curve then displays log-periodic oscillations with a tem¬ 
perature period set by the Efimov state energy spacing 
of exp(27r/so) « 515, where sq = 1.00624, and with a 
phase given by The relative peak-to-peak ampli¬ 

tude is 7% for ^^^Gs. As seen in Fig. 4, such oscillations 
cannot be resolved in the experimental data because of 
limited signal-to-noise and the limited range of tempera¬ 
ture. The predicted contrast of these oscillations for ^Li 
is even smaller (^ 6 %). This is a general property of the 
system of three identical bosons due to the smallness of 
|sii| [15]. 


IV. CONCLUSIONS 

In this article, we developed a general theoretical 
model for the coupled time dynamics of atom number 
and temperature of the 3D harmonically trapped unitary 
Bose gas in the non-degenerate regime. The theory takes 
full account of evaporative loss and the related cooling 
mechanism, as well as of the universal three-body loss 
and heating. It is furthermore extended to the special 
case of 2D evaporation. We predict and experimentally 
verify the existence of a “magic” trap depth, where the 
time derivative of temperature vanishes both in 3D and 
2D evaporation. 

We compare our model to two different set of experi¬ 
ments with lithium and cesium with vastly different mass 
and temperature ranges. The data illustrates the univer¬ 
sal scaling over 2 orders of magnitude in temper¬ 
ature, and we obtain an experimental value of the Efi¬ 
mov inelasticity parameter for the 47.8-G resonance in 
^^^Cs. The theory further enables an independent deter¬ 
mination of the trap depth. The agreement found here 



with standard methods shows that it can be used in more 
complex trap geometries (crossed dipole traps, or hybrid 
magnetic-optical traps) where the actual trap depth is 
often not easy to measure. 

In future work it would be highly interesting to probe 
the discrete symmetry of the unitary Bose gas by reveal¬ 
ing the 7% log-periodic modulation of the three-body loss 
coefficient expected over a factor 515 energy range. 
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Supplemental Material: Universal Loss Dynamics in a Unitary Bose Gas 


A. setup 

Our setup is a modified version of the one presented in [34]. The ^^^Cs atoms are trapped by means of three 
intersecting laser beams, and a variable magnetic field gradient in the vertical direction (partially) compensates 
gravity. An intrinsic advantage of the scheme is the perfect spin polarization in the lowest hyperfine ground state 
\F,mp) = |3, 3), because the dipole trap potential is too weak to hold atoms against gravity if they are in any 
other ground state. As we will see, the trap frequencies stay almost constant when reducing the trap depth, making 
evaporation very efficient [34] . 


Trap model 

The trap consists of three 1064-nm laser beams and an additional magnetic gradient field, see Fig. SI (a). All 
beams propagate in the horizontal plane. An elliptical light sheet beam (power Pls =520 mW, waists: wlsv = 
33.0 pm [vertical direction] and wush = 225 pm [horizontal direction]) creates the vertical confinement together with 
the magnetic gradient. Two round beams (Tr,i,r 2 =1.1/1.2W, waist of wr = 300 pm ) stabilize the horizontal 
confinement. The light sheet center is zq = 6 pm lower than the center of the trap formed by the round beams only. 
The potential along the vertical axis can therefore be written as 

V(z) = -ULse ”£sv-[/j^e +jz, (SI) 

where the Uls and Ur are the contributions from the light sheet and round beams, respectively. The tilt 7 has a 
gravitational and a magnetic contribution. 


J = mg- (S2) 

where m is the atomic mass, g is the gravitational acceleration, B' = dzB^ is the magnetic held gradient along the 
0 -axis, and p = 0.75pe is the atom’s magnetic moment in the |3,3) state, with pB being Bohr’s magnetic moment. 
Thus, a gradient of Bq — B'('y = 0) = 31.3 G/cm is needed for magnetically levitating the cloud. An example 
potential shape is given in Fig. SI (b). 



FIG. SI: (Color online) (a) Schematic drawing of the hybrid trap. It consists of three intersecting lasers, and a magnetic field 
gradient in the 0 direction (vertical) created by a pair of coils. The light sheet beam (blue) confines dominantly along the 
vertical direction. The additional round beams (red) stabilize the horizontal confinement, (b) Trap shape along the vertical 
( 2 —) direction. It is composed of two Gaussians and a linear contribution from the tilt, see Eq. (SI). Relative dimensions are to 
scale. Also indicated are the contributions of the tilt only (dashed line) and the round beams (dotted line), (c) Trap frequency 
measurements (dots) and fits (line) as a function of the tilt of the trap. The frequencies are normalized with respect to the 
zero-tilt frequency = 2n x 140 Hz. We normalize the tilt from our knowledge of the critical 7 c, where the trap opens, see 
text. 
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Trap frequency calibration 


When we intentionally change the trap depth, we also change the trap frequencies, mostly affecting the vertical 
direction. The data was taken during two different measurement campaigns in 2012 and 2013. Therefore, the trap 
had to be recalibrated for each of this campaigns, and the data is presented in a normalized way. 

We measure the oscillation frequency along the z axis as a function of the tilt, see Fig. SI (c). This is established 
by inducing sloshing oscillations to a small, weakly-interacting Bose-Einstein condensate (BEC), and performing 
time-of-flight measurements of its position after a variable hold time. 

We fit the measured 0 -axis frequencies to a numerical model of the trap potential Eq. SI. In the model, we plug 
the aspect ratio of the trap depth contributions from the three beams i. 




2 aPj 

TT WhWv ’ 


(S3) 


where a/ks = 2.589 x 10“^^K.cm^.mW“^ is the atomic polarizability at 1064nm [35], Pi is the power in beam i, and 
the Wh/v are the waists in the horizontal or vertical direction. We are left with two fit parameters: The frequency 
at zero tilt U!z(j = 0) = oJz^o, and the critical tilt jc = l{B' = B'^) where the trap opens (local minimum in V{z) 
disappears) and oj^ goes to zero by construction. We find oj^^o = 27r x 139(1) Hz and 27r x 140(1) Hz for 2012 and 2013, 
and P' = —0.3(4) G/cm and —4.3(7) G/cm for 2012 and 2013. With this calibration, we introduce the normalized 
tilt 7 / 7 c- The values for P' coincide well with the gradient values observed when increasing the tilt until a small 
(< 5000 atoms) weakly-interacting BEC drops out of the trap. 

The kink in the trap depth theory curve (Fig. S2) near 7 / 7 C = 0.02 (shaded area in Fig. S2(a)) corresponds to a 
situation depicted in Fig. SI (b), where the contribution of the large-waist horizontal beams on the trap depth vanishes. 
The blue-shaded region of the horizontal beams’ contribution extends over the small region from 7 / 7 C = —0.02 to 
0.02. Therefore, small experimental uncertainties on the applied magnetic gradient, or additional trap imperfections 
can explain the fact that we do not find this sudden rise in U. Other than that, we see a remarkable correspondence 
between theory and experiment. 

We note that the data can also be well described by the analytical model of a single gaussian potential (Pr = 0) 
with tilt, as presented in [34]. Because of the large mismatch between wr and wlSj the presence of the round beams 
mainly affects the horizontal trapping. The critical gradient we find is only 2% larger than the single-gaussian value 
V^7 Wls/2PlSv [34]. Furthermore, the horizontal trapping frequencies (ujx,ujy) « 27r x (13,30) Hz, measured with a 
similar method for each dataset, remain constant. 


Imaging system calibration 

The high-resolution imaging system is similar to the one presented in [14] . It is well calibrated using the equation 
of state of a weakly-interacting 2D Bose gas for the absorption-coefhcient-to-atomic-density conversion (in good 
accordance with the method of classical 2D gas atomic shot noise [36]). The imaging magnification is obtained from 
performing Bragg spectroscopy on a 3D BEG, using the variable retroreflection of the 1064-nm round beams. 


^®^Cs sample preparation 

We prepare the ^^^Cs samples in the trap described before. In brief, after magneto-optical trapping and degenerate 
Raman sideband cooling we obtain magnetically levitated (7 = 0) samples of 10® ^®®Cs atoms at 1 |xK [34]. We can 
cool the samples further by evaporative cooling. In order to achieve this, we adjust the trap depth U{'y) by changing 
the tilt 7 of the potential (SI). Thus, the samples can be evaporatively cooled all the way to quantum degeneracy in 
« 2 s [34] at 20.8 G, yielding a scattering length of 200 Og, with ug being the Bohr radius [29]. 

We prepare our samples by stopping the evaporation at a given tilt. We then ramp the tilt adiabatically to the 
desired value. Finally, at a time to, we jump the field to the Feshbach resonance at 47.8 G [26] in typically < 1ms 
and wait for at least ‘I'k jizix in order for the samples to reach dynamical equilibrium. We are therefore able to prepare 
samples of variable initial parameters: Atom number Wn, temperature Tin and relative trap depth rjin = U{'j)/ksTin, 
where fee is the Boltzmann constant. After a hold time t, we take an in-situ absorption image with a vertical imaging 
setup. 
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B. ^Li setup 

The ^Li data was taken using the apparatus described in [15]. This trap consists of a 1073-nm single-beam 
optical dipole trap providing adjustable radial confinement, and an additional magnetic field curvature providing 
essentially infinitely deep harmonic axial confinement along the beam axis. After loading into this trap, the gas is 
evaporated by lowering the radial trap depth at a magnetic field of 720 G, where the 2-body scattering length is 
200 oq. The evaporation in this hybrid trap is then effectively 2D. After the temperature and atom number of the 
gas have stabilized, the radial trap is adiabatically recompressed by about a factor of two. Since the axial magnetic 
confinement is practically unchanged, this recompression causes the temperature of the gas to increase with lot , 
while the trap depth increases as uP'. Consequently, by varying the amount of recompression, we can vary rjia- After 
this recompression, the magnetic field is ramped to the Feshbach resonance field of 737.8 G in 100-500 ms and N{t) 
and T{t) are measured with in situ resonant absorption imaging perpendicular to the long axis of the cloud. The trap 
shape can be described by Eq. (SI) with Dls = 0, wr = 37(1) pm, and replacing z by p, the radial coordinate. I/r is 
the dipole trap potential with power Ptrap and a is the polarizability of the ^Li atoms at 1073 nm. We can neglect 
the tilt 7 because of the small mass of ^Li. 


C. Time scale order 

In order for the theory to be valid, we make sure the timescale order is not violated: 


where we have the three-body loss time constant (cf. Eq. (19)) 

5 


^3B 


the evaporation time constant (cf. Eq. (25)) 


Tev = 772 1 e 






the two-body scattering time constant 


'^2B — n,Qcr(7U, 


and the trapping time constant 


Arap — Wslow) 


where Wsiow is the slowest trapping frequency (along z in the ^^^Cs case). 


(54) 

(55) 

(5 6 ) 

(57) 

(5 8 ) 


D. Fits to the model 

For each data set, we have decay data for N{t) and T{t). We fit both temperature and atom number individually 
with solutions to the coupled differential equation set of Eqs. (12) and (25). For both fits, we use a common three-body 
loss coefficient A 3 , and a common trap depth U. The fitting is done by minimizing the weighted sum axT + oi~^XN 
by varying both the weighing factor a and the fit parameters. The quadratic deviations are defined as xt,n = 

{o'T,N being the deviations of data and fit). This method also accounts for the different amount of relative signal-to- 
noise ratio of both data sets. 


E. Trap depth 

As an independent test of the theory fits, we compare the fitted trap depth U to its independently known counterpart 
from experimental parameters. In Fig. S2(a), we plot the ^^^Gs results as a function of the relative trap tilt x/lc- 
We also plot the theoretical value for U{x/xc) as a solid line. Except near zero tilt, we find excellent agreement 
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Relative tilt yly^ Power Pj^ap in W 

FIG. S2: (Color online) (a) Relative trap depth results from the fits to our data (dots), and theoretical model (Eqs. (SI) 

and (S3), line). The trap depth is normalized with respect to the trap depth Uls ~ 11 [rK given by the gaussian light sheet 
only. The shaded area corresponds to the region where the horizontal confinement beams significantly contribute to the trap 
depth, (b) Absolute ^Li trap depth results from the fits to our data (dots). The solid line indicatess theoretical knowledge of 
our trap (Eqs. (SI) and (S3), line), with wr = 38(1) qm. The shaded area accounts for the combined uncertainty of wr and 

Ttrap. 


of the fitted values with the values known from experimental parameters. For ^Li, see Fig. S2(b), we find excellent 
agreement with our theoretical knowledge of the trap depth, which is given by the dipole laser waist wr, power Ptrap 
and the atom’s polarizability. It is indicated by the shaded area in Fig. S2(b). Therefore, we can infer the dipole trap 
laser’s waist in an independent fashion. From the fit to our measured trap depths (solid line) in Fig. S2(b) we obtain 
Wr = 38(1) pm. This value coincides with independent measurements of wr from fitting the trap frequencies as a 
function of Ptrap- These results emphasize the validity of the theory model (Eqs. (12) and (25)). 







